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Abstract 

We study the influence of a point defect on the profile of a growing surface in the 
single-step growth model. We employ the mapping to the asymmetric exclusion 
model with blockage, and using Bethe-Ansatz eigenfunctions as a starting approxi- 
mation we are able to solve this problem analytically in two-particle sector. The dip 
caused by the defect is computed. A simple renormalization group- like argument 
enables to study scaling of the dip with increasing length of the sample L; the RG 
mapping is calculated approximately using the analytical results for small samples. 
For a horizontal surface we found that the surface is only weakly pinned at the 
inhomogeneity; the dip scales as a power law with 7 = 0.58496. The value of the 
exponent agrees with direct numerical simulations of the inhomogeneous single-step 
growth model. In the case of tilted surfaces we observe a phase transition between 
weak and strong pinning and the exponent in the weak pinning regime depends on 
the tilt. 
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1 Introduction 



There has been much activity in the field of surface growth in the last period 
[1,2]. Some problems are now rather well understood, nevertheless, many ques- 
tions still remain open or not completely clarified. One of them is the effect 
of defects on growth. In particular, we have in mind a defect which persists 
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during the growth and makes the growth inhomogeneous. This problem is 
relevant for various physical situations. For example, 1+1 dimensional epitax- 
ial growth (one-dimensional substrate and one dimension in the direction of 
growth) of steps on vicinal surfaces when there is a line defect on the terrace in 
the growth direction. Growth becomes inhomogeneous, the growth velocity at 
the defect can be lower (the moving step edge is pinned to the defect) or higher 
(the defect causes an excess of incorporated particles) than averaged velocity 
of the homogeneous step edge. An impurity floating on the growth front can 
have similar effect. In the case of 2-1-1 dimensional growth on a substrate with 
a defect (for example a dislocation) the defect is replicated in the grown mate- 
rial and makes the growth inhomogeneous. Still another situation is growth of 
a two-component material forming two domains separated by a domain wall; 
growth rate at the domain wall is different from growth rate within domains 
[3]. Epitaxial growth in these situations was so far httle studied theoretically, 
although it is of great practical importance. Nevertheless, similar problems 
were investigated in different contexts [4-9] . 

From the statistical-mechanical point of view the central problem is the ex- 
istence and the character of pinning - depinning transition. The lack of de- 
position may result in a dip of depth d located at the defect. The growing 
surface may be pinned, if d scales linearly with the sample size L, or it may be 
depinned, if d remains finite when the sample size goes to infinity. However, 
the dip may scale as a power oi d ^ U' . To distinguish between differ- 
ent situations wc shall call the case with 7=1 strong pinning, and the case 
< 7 < 1 weak pinning. The situation when 7 = corresponds either to de- 
pinning, if lim^^oo d = const < 00, or to various types of logarithmic pinning, 
like d ~ (logL)'' with positive exponent 1/ . 

Several authors investigated the problem of pinning of growing interfaces 
through various formulations. Wolf and Tang [4] studied inhomogeneous growth 
analytically using the growth equation proposed by Kardar, Parisi and Zhang 
(KPZ) [10]. They found that depending on the sign of the coefficient A in front 
of the non-linear term in the KPZ equation there is different behaviour of the 
surface profile with the size L of the sample for the lack and for the excess of 
deposition at the defect. For A > 0, extra deposition at the defect site leads to 
a pile of amplitude d increasing linearly with L, whereas the lack of deposition 
produce a groove of depth d, increasing logarithmically, d ~ log L. For A < 0, 
linear and logarithmic dependence are interchanged. 

Kandel and Mukamel [5] performed numerical simulations of 1-1-1 dimensional 
polynuclear growth model with a defect. They observed a defect-induced phase 

transition when changing the strength of the defect, manifested by the change 
of the dependence of the dip, d, on the system size. They found that d scales 
as d ^ L^, with 7 = 1 in the strong-defect phase, while 7 < 1 for weak defect 
phase. In further investigation, a mean-field description of the phase transition 
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was used [8] and the mean-field value of the exponent 7 in the weak defect 
phase was found to be 7 = 0. 

Tang and Lyuksyutov [6] investigated the problem in mathematically equiva- 
lent formulation as localization of a directed polymer in disordered medium. 
They found that there is no depinning transition in 1+1 dimensions for a line 
defect[^ and in 2+1 dimensions for a plane defect, while depinning occurs for 
a line defect in 2+1 dimensions. This finding is not in contradiction with the 
results of [5], because in [6] the short-time regime L 3> t'' is observed, where 
the weak pinning cannot be distinguished from strong pinning. 

The problem of growth in the single-step model is equivalent to diffusion in the 
asymmetric exclusion model [7], where the defect corresponds to a blockage. 
Such a situation was already studied analytically and some results on the 
nature of the impurity-induced phase transition are available. They deal with 
the presence of a strongly pinned phase. The matrix technique developed by 
Derrida et al. [7] for studying the asymmetric exclusion model with open 
boundary conditions was used in the case of mobile impurity in the medium, 
observing a phase transition when changing the strength of the impurity [11]. 

Using Bethe Ansatz it was also argued in favour of the existence of the phase 
transition induced by a point defect [12]. The pinning-depinning transition 
corresponds to the phase transition characterized by condensation of particles. 
The exact analysis of small systems {L < 4) accompanied by the use of Fade 
approximants [13] gives some rigorous bounds on the region of strongly pinned 
phase. However, the possibility of the phase transition to the weak pinning 
regime remained open. 

To our knowledge, the case of weak pinning was not explicitly considered 
in the existing literature and none of available calculations yields reasonable 
estimate of the exponent 7 in the weakly pinned phase, even though the matrix 
technique should in principle be able to give exact solution and the generalized 
Bethe Ansatz for complicated boundary conditions seems promising [14]. At 
present it is not clear under which conditions the weakly pinned phase with 
< 7 < 1 exists and what is the value of 7. 

In this paper we show the existence of weak pinning and present the calcu- 
lations of the exponent 7 in 1+1 dimensional single-step model of inhomoge- 
neous growth. We use an approach combining Bethe Ansatz and renormal- 
ization-group argument. We also present results of numerical simulations con- 
firming the predicted exponent. 



^ In the context of directed polymers the term line defect corresponds to what 
we shall call point defect in the growth model or asymmetric exclusion model 
terminology. 
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The article is organized as follows. In section 2 we describe the growth model 
considered, and give the definition of the dip for the surface of an arbitrary 
orientation. In section 3 the one particle solution is described explicitly. The 
latter result is used in section 4 for construction of two-particle states. Section 
5 contains our main results, both analytical and numerical, and the last, sixth 
section is devoted to discussion of the results. 



2 Growth model with a local defect 



In numerical simulations of surface growth, solid-on-solid (SOS) models are 
often employed because the surface can be simply described by a single height 
function h{x) of the substrate coordinate x. In most of numerical studies of 
inhomogeneous growth, so-called restricted solid-on-solid (RSOS) models with 
the additional constraint on the height difference of the neighbouring sites 
were used. The direct consequence of the constraint is that the number of 
possible configurations is finite. The RSOS models are rather well understood 
in the case of homogeneous growth [1]. The simplest one of these models is the 
single-step model with the height difference being only -|-1 or — 1. This class of 
models is particularly suitable for numerical as well as analytical approaches. 
In 1-1-1 dimensions an exact solution is known, using the Bethe Ansatz [15]. 

We consider an inhomogeneous 1+1 dimensional single-step growth model. 
The state of the system is described by an array h{x),x = 0,1,2,. of integers, 
which represent the height of the surface at the site x with the constraint 
\h{x + 1) — h{x)\ = 1 which enables representation of the state of the surface 
using Ising spin variables a{x) defined by a{x) — h{x) — h{x — 1). We impose 
cyclic boundary conditions. In terms of the spin variables they read (t(0) = 
cr{L), while for the height variables they fix the difference h{0) — h{L), i. e. 
the global tilt, constant during the growth. We consider only pure growth 
(no evaporation). Growth is possible only at a growth site x which obeys the 
condition a{x) — —1, a{x -|- 1) = 1. The deposition of a particle results in the 
change of the configuration: h{x) — > h{x) + 2 i. e. a{x) — >• —a{x), a{x -t- 1) — > 
—(t{x + 1). The inhomogeneity is introduced by the change of the rate for 
deposition on one fixed site. We suppose the deposition rate at a growth site 
X ^ xq to be equal to 1, and the deposition rate at a growth site on the defect 
site a; = Xo to be 1 — q; . The variable a denotes the strength of the defect. 
The value a — 1 corresponds to total inhibition of growth at the defect, i.e. 
to maximum possible pinning; due to the single-step constraint the surface 
profile evolves to a triangular shape in which it remains frozen. 

Let h{x,t) be the surface height at time t. To measure the pinning of the 
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surface with an arbitrary tilt {h{L,t) — h{0,t))/L, we introduce the quantity 

AM-) = ^ g (m-, t) - MO, t) - ,MM^MM) j , (1) 

which we call averaged height difference. It is the height averaged over (long) 
time interval T with the steady increment of the height and the tilt of the 
surface subtracted. If there were no defect, the averaged height difference 
would be zero. The defect at xq introduces an inhomogeneity into the system 
and Ah{x) may depend on position. We define the dip at the defect site as 

^^tEAMx) -A^xo). (2) 

^ x=l 

The surface is pinned, if d/L remains finite in the thermodynamic limit, L 
oo. It seems natural to say that the surface is depinned ii d/L — > 0. In this 
sense, Xuhl^^ is the order parameter of the pinning transition. However, 
d/L may converge to zero in the thermodynamic limit even if o? ~ L''' with 
< 7 < 1. The dip blows up, but more slowly than the length of the sample, 
therefore we call this situation weak pinning with exponent 7. 

The dynamics of the system has the following conservation law. The sum of 
all spin variables, S — Y.x '^(-) does not change during the growth. S should 
be understood as the average orientation of the surface, S — h{L) — h{0), 
fixed by the boundary conditions. The single-step model is mapped to the 
asymmetric exclusion model if we consider down spins as particles and up 
spins as unoccupied sites. In the language of asymmetric exclusion model, the 
corresponding conserved quantity is the number of particles n — {L — S)/2. 
The particle density p — n/L is 1/2 for horizontal and p ^ 1/2 for tilted 
surface. 

The situation corresponds to directed diffusion of n hard core particles on a 
line. At each time step, each particle can remain in the original position or 
jump to its right next neighbour site, provided it is unoccupied. The jumping 
rate from the site to + 1 is 1 — a, while the jumping rate from x io x + 1 
for X ^ XqIS 1. 

The conservation law facilitates the exact solution of the system. The Bethe 
Ansatz solution [15] starts with writing the master equation for the probabil- 
ities of all possible configurations. The set of all configurations can be divided 
into sectors according to number of particles. For each number of particles n, 
the master equation 

p{xi, X2,...,Xn) = Tn{x}{x'}P{x'i, ...,x'^) (3) 

{x'} 
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is solved independently. Here, the variables xi, ...,Xn describe positions of the 
particles and p{xi, ...,Xn) is probability of finding the particles in the spec- 
ified spatial configuration. The spectrum of the n-particle transition matrix 
Tji contains all information on the dynamics of the system. The solution for 
arbitrary n is still not known. We shall present exact solution for n = 1,2. 



3 One-peirticle problem 

In the case of a single diffusing particle the solution is quite easy. For n — 1 
the master equation has a very simple form 

p{x) = p{x - 1) - p{x) , for Xy^xo,Xy^xo + l , 

p{xo) = p{xo - 1) - (1 - a)p{xo) , (4) 
p{xo + 1) = (1 - a)p{xo) - p{xo + 1) . 

We can solve it by finding the (right) eigenvectors of the transfer matrix Ti 
in the one-particle sector. The eigenvectors are parametrized by a complex 
number C. As we shall see, also the left eigenvectors are useful, so we determine 
both of them. 

We write 

y 

for the right eigenvectors and 

Y,T^c{x)Ti^y = X(7i'^{y) (6) 



for the left ones. The left and right eigenvectors for different are mutually 
orthogonal. 

We look for the solution in the form 7r^(x) = A(^^ for x < xq or x > xq, where 
A may be different in the regions x < xq and x > xq. We find easily 

7r^(x) = , for X < , 

M^o) = e°T^ , (7) 
M^) = C'^iE^ , ioTx> xo , 
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and similarly for the left eigenvectors 



TT^ (x) = C ^ , for X < Xo , 

7rf(^) = C-^tE^ , for X > Xo . 



(8) 



The corresponding eigenvalue is = — 1. The stationary solution Ps{x) of 
the master equation is the eigenvector with zero eigenvalue, i.e. with C — 1, 
normalized so that J2xPs{x) = 1. The stationary solution gives immediately 
the stationary height profile. For a — 0, the value of Ah{x) is zero for all x, 
whereas for a > we have 



^H^) = L(L(l-a)+a) , ^01 X < Xq 
^^(^) = L(T(ta)la) , for X > Xo 



(9) 



Hence the dip at the defect is 



Because there is only one particle, it is obvious that the surface is not pinned; 
the value of the dip for all a and L is at most 1. So, in order to actually 
observe the pinning, we should have non-zero density of particles in the ther- 
modynamic limit L — i> cxD, n — > oo. This requires a substantially different 
approach. 



4 Lattice constant doubling 



We apply a scheme similar to the real space renormalization group technique. 
Suppose we have even number n of particles in the sample of length L, which is 
also even. The time evolution is governed by the transfer matrix T„(q;). (From 
now on we shall write explicitly the dependence of the transfer matrix on a). 
The presence of the defect of strength a results in a dip d{n, L, a). Taking each 
pair of neighbour particles as a single object, we have system of n/2 effective 
particles. Their interactions are rather complex, but the fact, that the parti- 
cles are hard-core, i.e. they can not skip one another, is preserved even for 
the effective particles. The time evolution is no more described by a Markov 
process governed by some transfer matrix, but by a more general evolution 
operator 7^/2 («)• However, we suppose that we can approximate the dynam- 
ics by a Markov process and map the system of n/2 effective particles onto 
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the original system, but with half number of particles and different effective 
strength of the defect a'2,: 



(11) 



The lattice constant of the new system is twice as large, because pairs of par- 
ticles occupy two original sites, so the length of the rescaled sample is half 
of the original one. Also the height is rescaled by a factor 2. The procedure 
roughly corresponds to replacing four atoms deposited independently on the 
surface by a single "composite" atom made of four original atoms glued to- 
gether. The dependence of rescaled strength of the defect 0:2 on the original 
a can be computed by equating the dip in the original and rescaled system 



This equation may be understood also as the requirement of the equality of 
deposited mass in the original system and the rescaled system. 

A justification for the doubling procedure may be provided a posteriori by the 
comparison of the whole surface profiles in the original and rescaled systems. If 
the deviation is small even far from the location of the defect, we can suppose 
that the pairs of particles actually behave like the original particles, but in 
rescaled system. 

The simplest way for obtaining an approximation for the function 02(0;) is to 
compute the average height profile for n — 2 and the corresponding dip. From 
(10) and (12) it follows that 



Fortunately, the n = 2 case is soluble exactly. If a = 0, the solution is given 
simply by Bethe Ansatz. The many-particle eigenvector is the linear combi- 
nation of products of one-particle eigenvectors. However, doing the same for 
a does not work, because the system lacks translational symmetry. But 
the Bethe Ansatz eigenfunctions pb may serve as a good starting point for 
further calculations, because they must obey the equation 



in all cases when both particles are far from the defect. Moreover, if we con- 
struct the eigenfunction from the exact one-particle eigenvectors (7), which are 
calculated for given a 7^ 0, then the equation (14) will hold even if one of the 




(12) 




(13) 



T2{a)pB = XpB 



(14) 
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particles is at the defect. Only such terms when both particles are at (or near) 
the defect will disagree. But we may treat the difference as a localized pertur- 
bation and solve exactly the equation for the resolvent G{z) = {z — T2)~^. 

So, let us start by exphcitly writing the eigenvalue equation (14). We define 
the hopping rate from the site x 

Wx^l -aS^x^ . (15) 

Then, the two-particle eigenvalue equation is 

A p{x, y) = p{x-l,y)+ Wy-i p{x, y - I) 

-{wr, + Wy)p{x,y) , (16) 
A p{x, X + 1) = Wx-ip{x — 1, a; -h 1) — Wx+ip{x, x + 1) . 

Prom the cyclic boundary conditions it follows that p{x, y + L) — p{y, x). 
We construct the Bethe eigenfunctions 

PbcU^: y) = 7r^(x)7ra,(y) + ABTloj{x)Tl(:{y) , (17) 

using the one-particle functions (7). The constant 

Ab = (18) 

is set so that (16) are satisfied far from the defect [15]. For C, —uj — Iwe. take 
Ab — 0. The corresponding A is 

A = Ascc. = C"' + w"' - 2 . (19) 

The values of C, and uo are set by the cyclic boundary conditions, resulting in 
a pair of algebraic equations 

(-L ^ C-i 1-qC 

^ ^-1 l-a ' ^20) 

, ,L io—l \—au> 

Q—1 1—a ' 

which can be solved numerically by standard methods. Together with the 
solution ( = u = 1 we have L(L — l)/2 pairs (C,'^)- Each of them determines 
one of right eigenvectors of some matrix T2B, which differs from T2 only in 
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several matrix elements, so we can take T2B as a first approximation to T2 and 
treat the difference T2—T2B as a perturbation. Let P be matrix whose columns 
are the Bethe eigenfunctions PBCui{x,y) (pairs {x,y) denote rows, pairs (C,'^) 
denote columns). We can write 

T2B = P-'AP (21) 



where A is the diagonal matrix with numbers A^i^t^ on the diagonal. It is quite 
easy to find the inverse if o; = 0. As Bethe functions pb are exact right 
eigenvectors of T2, the Bethe functions p]^^^ constructed by the same procedure 
from the left one-particle eigenvectors tt^ are exact left eigenvectors of T2. 
Because left and right eigenvectors with different eigenvalues are orthogonal, 
the matrix composed of properly normalized left eigenvectors is the inverse of 
P. When a 0, the situation is different, but in fact we need to make only a 
tiny modification. 

It can be checked directly (actually we performed the check using Maple V) 
that the vectors 

PB(x,y) ^^I^(x)7rl(y)+ABnf(x)7r^(y) , 

for {x,y) {xq,Xo + 1) , (22) 
pI{xo,xo + 1) = + AbC^'^uj-'^"-') 



are orthogonal to the vectors pb and so, after proper normalization, they 
represent rows of the matrix P~^. 

As we already said, vectors ps are considered as nearly exact eigenvectors of 
T2. What we exactly mean is that 

T2P = PA + A , (23) 



where the matrix A has only two non-zero rows, corresponding to row indices 
a = {xq — 1, xq) and b — {xq, xq + I). (The column indices are the pairs {(, lu)). 
We have explicitly 

.,-1 _ T ^^Ai^ 





= «(c-^ 




= -a{l 


Aa(i,i) 


a 

~ 1-a ' 


A6(i,i) 


a 

~ 1-a 



C-1 (l-aC)(l-Q;a;) ' 
c-1 (l-aC)(l-au;) ' ^24) 
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The solution is obtained, as soon as we compute the resolvent 
G{z) = {z — T2)^^. If we define the "unperturbed" resolvent G°{z) = {z — 
P~^AP)~^, using elementary matrix algebra we obtain equation 

G{z)^G"{z) + G{z)A{z-A)-^p-^ , (25) 

which we call Dyson equation, taking an analogy from the quantum mechanics. 

Let n be projector to the two-dimensional subspace corresponding to indices 
a and b and n = 1 — H the complementary projector to it. The Dyson equation 
can be solved within that subspace, because 11 A = A. We obtain 

G{z)U = G"{z)U{l - A{z - A)-ip-in)-i 

G{z)U = G°(^)n + G°{z)U{l - A{z - A)-ip-in)-i x (26) 

xA(z- A)-^p-^n . 

Computation of the matrix inverses is very simple, because they involve only 
matrices 2x2. 

The stationary state is obtained from the resolvent in the ^ — > limit 

Ps{x, y) ^\im zG^^,y)a(z) . (27) 

The column index of the resolvent may be chosen arbitrarily. It corresponds 
to the fact that stationary state does not depend on initial conditions. 



5 Results 

We calculated explicitly the two-particle (n = 2) stationary state Ps{x,y), 
defined on the triangle I < x < y < L for L ranging from 4 to 32. The figures 
1 and 2 show the stationary state for L = 12 and a — 0.1, 0.9. In both cases, 
the defect is located &t xq — 2. We can see that for small a the effect of the 
defect is weaker, but more delocalized in the two-dimensional configuration 
space {x,y). The probability is distributed nearly homogeneously, but there 
is a band of higher probability composed of points which have any of the 
coordinates x, y equal to xq. Large a leads to localization of the eigenfunction 
around the point {xq — 1, xq); in the latter case a single peak develops which 
bears nearly all the probability. 

These results were further used to calculate the stationary height profile differ- 
ence Ah{x) and dependence a2{a) of the renormalized strength of the defect 
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on the bare value of a. The figure 3 shows the profiles for L = 4 and different 
a's. This situation corresponds to horizontal surface (p = 1/2). In the same 
Fig. 3 we can also compare the exact n = 2 (two-particle) profiles with the 
effective one-particle ones, calculated for n = 1, L = 2 and defect strength 
0:2(0;). We can see that the agreement is rather good in the whole range of 
a, which justifies our doubling procedure. The figure 4 shows the same for 
L = 12. These height profiles are tilted (p < 1/2) which leads to the asymme- 
try of the dip. Even in this case we can see good agreement between the exact 
two-particle profile and the effective one-particle one. However, for small a 




Fig. 1. Stationary state of the two-particle system for L = 12 and a = 0.1. 




Fig. 2. Stationary state of the two-particle system for L = 12 and a = 0.9. 
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and near the defect, the agreement is worse than for horizontal surface. 



The key result of our computation is the function a2{a). In the figure 5, the 
dependence of a2{a)/a on a is shown for both horizontal (p = 1/2) and tilted 
(p < 1/2) surfaces. 

The physics of pinning is contained in fixed points of the map a a2{<y). 
There are surely at least two fixed points, namely and 1, corresponding to 
weak pinning (or depinning) and strong pinning regime, respectively. Their 
properties may be determined from the quantities 



Sn — limn,_>n ~ , 

^ ' (28) 
si = lim^,^! . 

The inequality Sj < 1(> 1) means that the fixed point a — i is stable (unsta- 
ble). 

Let us comment on the weak pinning regime first. The value of Sq determines 




Fig. 3. Stationary height difference for the two-particle system for L = 4 and 
a = 0.1,0.3,0.5,0.7,0.9 (solid lines) together with effective profile of corresponding 
one-particle system (dashed lines). The numbers close to each full line denote the 
corresponding a. 
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the scaling of the dip in the weak pinning regime for o; — > 0. The exponent is 



as can be seen from equations (10) and (12). 

The values we observe in the Fig. 5 are 1 > sq > 0.5 in all cases, corresponding 
to weak pinning, 1 > 7 > 0, while Sq < 0.5 would be signal of depinning (or 
logarithmic pinning), 7 = 0. Thus, the weak pinning regime is always present. 

The most important result is the value of the exponent for horizontal profile. 
We have observed an interesting result, that for p = 1/2 it is Sq = 3/4 up to 
the precision of our numerical procedure (10~^). So, we conjecture, that the 
exact value of the exponent for horizontal surface and a — > is 

7 = = 0.58496... . (30) 

' ln2 ^ ^ 



This value may be compared directly with numerical simulations. 




Fig. 4. Stationary height difference, for the two-particle system for L = 12 and 
a = 0.1,0.3,0.5,0.7,0.9 (soHd Unes) together with effective profile of corresponding 
one-particle system (dashed lines). The numbers close to each full line denote the 
corresponding a. 
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The dependence of sq and the corresponding exponent 7 on the particle density 
p = n/L is shown in Fig. 6. Wc can see that the exponent 7 decreases smoothly 
when decreasing the density and seems to approach for p — > 0. It reflects the 
fact, that for zero density there is only finite number of particles on infinite 
sample and the dip should be also finite. The point at p— 1/2 does not fit to 
this smooth dependence. The explanation lies most probably in the different 
symmetry group of the horizontal surface. (For tilted surfaces, the particle-hole 
symmetry is broken.) 

Let us now turn to the discussion of the phase transition from weak-pinning 
to strong-pinning regime. 

From Fig. 5 we can see, that ct = is stable fixed point because sq < 1 for 
all p. On the other hand, we found that the fixed point o; = 1 is unstable for 
p = 1/2 (si > 1), but stable for p < 1/2 (si < 1). 

The existence of a phase transition for some a = ac would be signalled by the 
presence of a third (unstable) fixed point. We may formulate it by saying that 
R{a) = 0:2(0;) — a has a root at < etc < 1 with positive derivative, -R(oc) = 0, 
Sc — R'{oic) > 0. The critical value ccc separates the weak- pinning (low a) and 




Fig. 5. The a-dependence of the ratio of the effective defect strength to original 
a for p = 1/2 (solid Une) and p = 1/3,1/4,1/5,1/6,1/8,1/10,1/12,1/16 (dashed 
Unes, from top to bottom). In the inset, detail of the dependence of difference a2 — a 
on the distance from the critical point, a — ac, for p = 1/3 (solid), p = 1/4 (dashed), 
p = 1/5 (dotted). 



15 



strong pinning (high a) regimes, corresponding to different scahng of the dip. 
For q; < we have 7 < 1, for o; > ccc we have 7 = 1. 

For p— 1/2 we did not find the third fixed point. Instead, the unstable fixed 
point at q; = 1 plays its role and we conclude that there is no phase transition 
for a < 1 and we are in the weak pinning regime for all a < 1. 

For p < 1/2 we find etc < 1, but very close to 1. The values of etc are 0.9371, 
0.9605, 0.97535 for p = 1/3, 1/4, 1/5, respectively. The dependence of R{a) on 
a — etc near the critical point is shown in the inset in Fig. 5. 

Finally, we compare some of our results with numerical simulations. We per- 
formed direct Monte Carlo simulations of growth in the inhomogeneous single- 
step model for horizontal surface, i.e., density p = 1/2 in asymmetric exclusion 
model. We measured the dip (2) as the function of the system size L and time 
t for several strengths of the defect a. When we start from the flat configura- 
tion, the dip increases with time and after some time tgat (growing with the 
system size) the stationary regime is reached. In this regime the dip fluctuates 
around the characteristic value dsat{L) ~ L"'. We concentrate on the weak 
defect regime, i.e., small a. In this case fluctuations are very strong and many 
independent runs are needed to obtain the value of the saturated dip. Results 
presented in Fig. 7 are averages over 10^ to 5 x 10^ independent runs. 

Wc observed that the slope, i.e. exponent 7 is slightly decreasing with de- 
creasing a, but it is deflnitely larger than 1/2. We calculated exponent 7 by 
fltting the last flve points in each curve. Values obtained for a = 0.05, a = 0.1 
and a = 0.2 are 0.59 ± 0.03, 0.61 ± 0.01 and 0.703 ± 0.004 respectively. The 
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Fig. 6. The slope so at a — (-I-) and the exponent 7 (x) against the particle 
density p = n/L. 
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value obtained for the lowest value a — 0.05 agrees well with the analytical 
prediction (30), so we conclude that the o; — > limit is well described by our 
analytical calculations. 



6 Conclusions and discussion 



We solved exactly in a special case the problem of growing interface in presence 
of a point defect. In the language of asymmetric exclusion model, we obtained 
the solution for cyclic boundary conditions in the two-particle sector. The so- 
lution consists in two steps. At first, using the Bethe-Ansatz eigenfunctions 
as a basis set, we found that the perturbation becomes localized in the con- 
figuration space. Then we solved the Dyson equation exactly by inversion of 
small matrices. We computed the stationary profile of the surface and deter- 
mined the dip at the defect. The two-particle result was then employed in 
a renormalization-group consideration, based on a lattice constant doubling 
transformation. We obtained the dependence of the effective defect strength 
a2 on the original, bare defect strength a. Fixed points of this mapping were 
identified. 

Our main result is the existence of the weak pinning phase and the value of 
the corresponding exponent. Let us discuss first the situation for a horizontal 




Fig. 7. Saturated dip, dgax, as a function of the system size, L, calculated in the 
inhomogcneous single-step model for three strengths of the defect a = 0.05, a = 0.1 
and a = 0.2 (from bottom to top). The dashed lines are fits to last five points, 
the slope of the full line corresponds to the value of 7 obtained by the analytical 
solution, eq. (30). 
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surface (p = 1/2). Wc have found just two fixed points: a = whicli is 
stable, and corresponds to weak pinning of the surface to the defect, d ^ 
with < 7 < 1, and a = 1 which is unstable. Hence, from our analytical 
calculations we obtained for a < 1 only the weak pinning phase. The strong 
pinning regime is limited to a single point a — 1. We calculated the exponent 
7 in the weak pinning phase to be 7 = 0.58496. 

We verified the existence of weak pinning by a direct Monte Carlo simulation 
of the inhomogeneous single-step growth model. For weak defect, a close to 0, 
we obtained 7 = 0.59 ±0.03 in good agreement with the analytical prediction. 

The Ref. [4] reports linear dependence of the dip on L, suggesting that 7 = 1. 
In fact, the authors of [4] state that d ~ sL, where s is the slope of the interface 
at the defect. They use the parameter s as a measure of the strength of the 
defect. We understand this choice as being dictated by the use of a continuous 
model (KPZ equation). However, we use different paramctrization, namely 
taking the strength of the defect a as a control parameter. For a fixed, the 
value of s may depend on L. A direct comparison of Ref. [4] with our results 
would require measurement of the quantity s, which was not performed in our 
simulations. 

In the case of tilted surfaces, we found that the fixed point a = 1 becomes 
stable, and in addition there is a third, unstable fixed point < ccc < 1 
which corresponds to the phase transition between strong and weak pinning. 
However, the value of etc in our calculations is very close to 1. 

The absence of phase transition at a < 1 for p = 1/2 is in contradiction 
with numerical simulations on the inhomogeneous polynuclear growth model 
[5] as well as with analytical results [12,13]. The most probable reason lies 
perhaps in the approximation used in calculation of the map «2(<^)- instead of 
working with macroscopic sample with generic length L and particle number 
n and compare it by the doubling procedure to the sample with length L/2 
with n/2 particles, very small sample, L — A was used. With such small 
sample, the influence of boundary conditions may be very strong, deforming 
the renormalization map 0:2 (ct)- Nevertheless, the agreement of the analytical 
results with simulations in the weak pinning phase shows that for weak defects 
the method we used works well. 

In order to improve the result, it would be necessary to solve the dynamics 
for a larger number of particles. However, the use of the Bethe-Ansatz eigen- 
functions for number of particles larger than two does not lead to locahzation 
of the perturbation and Dyson equation cannot be solved exactly by explicit 
matrix inversion. 

It would be of interest to investigate the existence of weak pinning in different 
inhomogeneous models, in particular one would like to know if the value of 
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the exponent 7 for the strength of the defect going to zero is universal, and 
eventually if a weak-pinning to depinning transition exists. 
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